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ABSTRACT 


The  molecular  cluster  model  approach  for  the  description  of  processes  at  surfaces  is  presented.  It 
involves  using  fmite  clusters  which  contain  surface  atoms  interacting  with  atomic  or  molecular 
adsorbates.  Accurate  ab  initio  wavefunctions  can  be  obtained  for  the  clusters  and  thus  makes  it 
possible  to  describe  the  bonding  in  terms  that  are  familiar  for  molecular  systems.  The 
adsorbate-substrate  interaction  can  be  anal>v.cd  by  means  of  new  techniques  developed  to 
characterize  the  bonding  as  covalent  or  ionic.  I'his  is  shown  for  F  on  Ag(l  1 1)  and  O  and  CO  on 
Cu(lOO)  systems.  The  measures  of  adsorbate  ionicity  arc  (a)  the  expectation  value  of  a  projection 
operator  giving  an  indication  of  the  number  of  electrons  associated  with  the  adsorbate,  (b)  the 
analysis  of  the  dipole  moment  curve  as  function  of  the  distance,  and  (c)  the  effect  of  an  c.xtemal 
electric  field  on  the  adsorbate  gcomctiy’.  Fhc  con.scqucnccs  of  ionic  adsorbates  on  metal  surfaces 
is  discussed  for  two  selected  topics:  adsorption  of  thiocyanate  on  Ag  electrodes  and  coadsorption 
of  CO  with  ionic  .species. 


.  Introduction 


The  concept  of  chemisorption  is  a  key  to  the  understanding  of 
the  catalytic  reactions  on  solid  surfaces.  Heterogenous  catalysis 
consists  of  elementary  reactions  occuring  at  the  surface  of  the 
catalyst.  In  these  reactions  chemical  bonds  are  formed  between 
the  surface  atoms  and  the  adsorbed  molecules.  These  interaction 
cause  rupture  of  chemical  bonds  within  the  adsorbing  molecule  and 
formation  of  new  bonds  between  the  resulting  fragments  and  the 
surface.  Thus,  the  first  step  for  the  understanding  of  catalysis 
is  understanding  of  chemisorption  phenomena. 

The  questions  that  are  asked  and  answered,  both  experimentally 
and  theoretically,  about  bonding  at  surfaces  are  very  similar  to 
those  that  are  asked  in  molecular  systems:  bonding  nature,  bond 
distances,  binding  energies,  vibrational  frequencies,  etc.  The 
theoretical  evaluation  of  these  energy-related  observables 
requires  the  computation  of  potential  energy  hypersurfaces  {PHS, 
often  called  potential  energy  surfaces  and  abbreviated  as  ?ES)  or 
of  their  sections  along  a  particular  direction  (potential  energy 
curves,  PEC)  for  the  interaction  of  the  adsorbed  molecule  and  the 
solid  surface.  However,  other  properties  are  often  easily  as 
interesting  as  the  total  energy  of  the  system;  for  example, 
dipole  moments  provide  useful  information  about  the  nature  of  the 
chemical  bond,  and  the  determination  of  dipole  moment  surfaces 
(curves)  is  required  for  the  evaluation  of  infrared  intensities. 

From  the  theoretical  point  of  view,  a  surface-adsorbate  system 
represents  a  very  complex  quantum-mechanical  many  particle  object 
which  requires  a  large  number  of  approximations  in  order  to  be 
treated.  Various  concepts  and  models  have  been  developed  for  the 


determination  of  electronic  wave  functions  describing  the 
interaction  between  the  adsorbate  and  the  substrate  [1]. 

The  cluster  model  approach 

In  the  cluster  model  approach,  a  finite  number  of  atoms  is 
chosen  to  model  a  local  site  on  the  crystal  surface  [2] .  In  this 
model  one  focusses  on  the  local  nature  of  interaction  between  the 
adsorbate  and  a  the  solid  substrate  which  involves  primarily  the 
nearest  neighbor  surface  atoms.  The  coupling  of  the  cluster  wave 
function  to  the  rest  of  the  substrate  (embedding)  would  have  to 
be  accounted  for  in  a  second  step.  Efforts  have  been  made  in  this 
direction  [3] .  Alternatively,  one  could  increase  the  size  of  the 
cluster  by  adding  additional  shells  of  substrate  atoms  until  the 
chemisorption  properties  become  independent  of  cluster  size.  When 
only  a  single  adsorbate  is  considered,  the  large  substrate 
cluster  represents  the  limit  of  zero  coverage,  \.e.  in  absence  of 
lateral  interactions  with  other  adsorbates.  Both  embedding  and 
cluster  convergence  are  formally  and  computationally  very 
difficult  probelms  so  that  often  isolated  clusters  of  small  or 
medium  size  have  been  studied  with  ab  initio  Hartree-Fock 
techniques  [4] ,  The  use  of  simplified  versions  of  the  Hartree- 
Fock  method  [5]  or  of  techniques  where  the  exchange  potential  is 
replaced  by  an  approximate  expression  depending  on  the  electron 
density  [6]  has  allowed  the  treatment  of  larger  systems  at  the 
cost  of  reducing  the  reliability  and  accuracy  of  the  results.  It 
must  be  noted  that  almost  all  theoretical  approaches  usually 
employed  in  quantum  chem.istry  have  been  applied  to  the  study  of 
cluster-adsorbate  systems. 

Since  in  the  cluster  model  approach  one  tries  to  simulate  a 
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solid  surface  by  means  of  a  finite  size  cluster,  one  usually  uses 
cluster  geometries  (and  eventually  electronic  states)  which  are 
appropriate  for  bulk  and  surface  geometries  of  the  substrate. 
These  clusters  have  little  in  common  with  the  real  geometries  and 
electron  distributions  of  free  metal  clusters  [7].  Indeed,  gas- 
phase  clusters,  often  highly  fluxional,  exhibit  a  very  cluster- 
dependent  chemical  reactivity  which  is  in  general  difficult  to 
rationalize  in  simple  terms  and  which  is  different  from  that  of 
bulk  materials  [8] .  Also  theoretical  computations  show  an 
oscillating  behaviour  of  the  chemisorption  properties  with 
increasing  cluster  size,  at  least  for  clusters  of  less  than  50 
atoms  [9] . 

The  big  advantage  of  the  cluster  model  approach  to  the  study 
of  chemisorption,  and  eventually  of  surface  chemical  reactions, 
is  that  the  interaction  can  be  analyzed  by  means  of  standard 
quantum  chemical  methods.  The  local  minima  of  the  total  energy 
Etot(Ri,  Rz  ,  R3  ,  ...)  as  function  of  the  nuclear  coordinates 
define  possible  equilibrium  geometries  of  the  cluster.  Usually, 
the  geometry  of  the  substrate  is  kept  fixed,  giving  a  static 
representation  of  the  surface;  however,  approaches  to  the 
chemisorption  of  atomic  or  molecular  species  in  which  the 
(partial)  relaxation  of  the  surface  has  been  taken  into  account 
have  been  reported  [10] . 

Cluster  wave  functions 

In  this  work  ab  initio  LCAO-Self  Consistent  Field  (SCF)  wave 
functions  are  used.  The  MO ' s  are  expanded  in  terms  of  Contracted 
Gaussian  Type  Orbitals  (GTO)  of  double  zeta  quality  or  better. 
Effective  Core  Potential  (ECP)  operators  [11]  are  employed  to 
describe  the  inner  shells  of  heavy  atoms  or  of  som.e  of  the 
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cluster  ruetal  atoms  surrounding  the  chemisorption  site.  Of 
course,  accurate  values  for  many  properties  cannot  be  obtained 
unless  correlation  effects  are  included  in  the  wave  function, 
e.g.  through  the  use  of  extended  configuration  interaction  (Cl). 
This  is  necessary  in  particular  for  the  exact  determination  of 
PES  (or  curves)  and  related  properties. 

Despite  the  limitations  due  to  the  cluster  size  and  trtfr'  usilllPi 
uncorrelated  SCF  wave  functions,  the  MO  cluster  model  can  provide 
valuable  information  about  surface  chemical  reactions.  In 
particular,  the  MO  cluster  model  is  ideally  suited  for  chemical 
interpretation  of  the  bond  at  a  surface  in  a  way  analogous  to 
that  used  for  molecular  systems.  Determination  of  SCF  PES ' s  is 
important  therefore  as  a  guide  in  identifying  mechanisms  which 
explain  the  surface  bonding  and  in  the  use  of  these  mechanisms  to 
correctly  interprets  the  significance  of  observable  properties 
and  dynamical  processes. 

Ionic  versus  covalent  adsorbates 

A  very  important  aspect  is  the  determination  of  the  ionic  or 
covalent  nature  of  the  adsorbate-surface  chemical  bond.  This 
question  is  considerable  fundamental  importance.  The  amount  of 
charge  transfer  from  the  surface  to  the  adsorbate  (or  viceversa) 
dramatically  influences  the  surface  reactivity.  We  mention  here  a 
few  selected  examples.  There  is  enough  evidence  that  the 
formation  of  strong  ionic  bonds  on  metal  surfaces  can  cause 
substantial  surface  reconstruction  [12];  also  the  energy  profile 
for  the  penetration  of  an  atomic  adsorbate  below  a  solid  surface 
can  be  strongly  dependent  on  the  ionic  character  of  the  bonding 
[13];  the  determination  of  the  charge  associated  to  chemisorbed 
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alkali  nietal  atoms  on  the  surface  of  semiconductor  materials  is 
of  basic  importance  in  order  to  understand  the  significance  of 
surface  metallization  and  of  the  metal-insulator  transition  [14] ; 
ionic  or  covalent  adsorbates  exhibit  completely  different 
behaviour  on  metal  electrodes  when  an  external  potential  is 
applied  [15];  finally,  the  fascinating  effect  of  promotion 
(inhibition)  in --catalysis  is  probably  connected  to  the  formation 
on  the  surface  of  ionic  species  [16] . 

The  identification  of  the  surface  chemical  bonding  as  covalent 
or  ionic  is  thus  of  fundamental  importance  in  surface  chemistry. 
An  analysis  of  the  adsorbate-substrate  binding  is  usually  based 
on  atomic  population  analysis.  These  are  often  misleading  because 
of  the  arbitrary  partitioning  of  electrons  between  different 
nuclear  centers. 

In  this  paper  we  present  new  techniques  developed  to 
characterize  the  chemical  bonding  at  surfaces  and  we  analyze  the 
consequences  of  ionic  surface  bonding  for  two  selected  topics. The 
paper  has  been  divided  into  two  parts.  In  Section  2  we  describe 
the  techniques  used  to  characterize  the  surface-adsorbate  bonding 
for  both  atomic  (F  and  O)  and  molecular  (CO)  adsorbates.  In  Sect. 

3  we  discuss  a  problem  relevant  to  electrochemistry,  the 
chemisorption  of  thiocyanate  ions  on  Ag  electrodes,  and  we 
present  a  new  mechanism  to  explain  the  effect  of  coadsorbates 
promoting,  or  inhibiting,  the  catalytic  activity  of  a  metallic 
substrate . 

2.  Analysis  of  the  surface  chemical  bonding 

To  determine  the  character  of  the  bond inn  of  an  atomic  or 
molecular  adsorbate  as  covalent  or  ionic  we  apply  three  measures. 
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No  single  measure  is  sufficient  to  unambiguously  characterize  the 
bonding,  but  the  three  measures  together  can  provide  a  reliable 
description  of  the  interaction.  The  analysis  is  applied  to  the 
following  problems:  a)  chemisorption  of  F  on  a  Ag^  cluster  model 
of  the  three-fold  symmetric  site  on  the  Ag  (111)  surface;  details 
of  the  cluster  geometry,  basis  set  and  ECP  used  can  be  found  in 
Ref.  17b)  Interaction  of  atomic  oxygen  with  d  Cus  ( 4 , 1  )••  cLuarih^ 
model  of  the  four-hollow  site  of  a  Cu(lOO)  surface  containing 
four  atoms  in  the  first  and  one  atom  in  the  second  layer;  the 
Cua -O  cluster  has  been  treated  at  all  electron  level  using  a 
[6s5p4d]  [18c]  and  a  [4s3p]  [18a]  contracted  GTO  basis  sets  for 
Cu  and  0  atom.s,  respectively,  c)  Interaction  of  CO  on  Cus  and 
Cui  4  cluster  models  of  the  Cu(lOO)  surface  (details  on  Ref.  11a). 
Projection  operator 

The  first  technique  makes  use  of  a  projection  operator  for 
some  or  all  the  occupied  orbitals  of  the  adsorbate  [19] .  The 
expectation  value  of  this  projection  operator  using  the 
wavefunction  of  the  whole  cluster-^-adsorbate  system  gives  a 
measure  of  the  extent  to  which  the  adsorbate  orbitals  are 
contained  in  the  complex.  Actually,  three  cases  are  possible: 
a)  the  expectation  value  Pp  of  the  projection  operator  P(p)  for  a 
given  orbital  p  is  close  to  2,  than  p  is  completely  occupied  in 
the  cluster  wavefunction  45;  b)  ?o  is  close  to  zero,  thus  p  is 
unoccupied  in  45;  c)  Pp  is  between  0  and  2,  indicating  that  a 
covalent  bonding  is  formed.  The  sum  of  ?o  over  all  orbitals  of 
physical  interest,  Np ,  provides  a  measure  of  the  number  of 
electrons  associated  with  the  adsorbate.  An  e.xample  of  this 
analysis  is  given  in  Table  1  for  F  on  Ag4  ,  0  on  Cu3  ,  and  CO  on 
Cuj  .  The  data  show  that,  for  distances  near  r.>  ,  the  ionicity  of 
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the  F  atom,  defined  as  Qp  =  Z  -  Mp ,  where  Z  is  the  nuclear 
charge,  is  very  close  to  -1.  F  is  close  to  be  a  perfect  anion. 
For  0  on  Cua ,  the  2s  and  2pa  orbitals  are  fully  occupied,  while 
P{2pm)  is  3.37  (Table  1).  In  this  case  Lt  is  necessary  to  take 
into  account  the  overlap  between  the  orbitals  of  the  isolated 
adsorbate  and  substrate  when  these  two  subunits  are  superposed  at 
r=re  .  This  overlap  amounts  to  0.41  and  therefore  the  2pn 
occupancy  is  intermediate  between  3.46  and  3.87.  The  oxygen 
ionicity  can  be  estimated  to  be  about  -1.6;  the  bonding  is 
covalent  with  large  ionic  character. 

The  analysis  for  the  last  case,  CO  on  Cus  ,  shows  that  the  five 
6  and  the  l.t  orbitals  of  CO  are  fully  occupied  in  Cus-CO.  The 
occupancy  of  the  CO  2-‘  MO,  1P(2-t*),  has  been  corrected  by  the 
overlap  of  the  two  subunits  and  represents  a  measure  of  the 
increase  in  the  ligand  2m*  occupation  due  to  metal  back-donation 
(Table  1).  P(2m)  depends  on  cluster  size.  The  projection  is 
consistent  with  the  CO  bond  being  covalent,  dominantly  in  m  space 
[11a]  . 

Dipole  moment  curves 

The  second  measure  of  bond  character  is  based  on  the  slope  of 
the  dipole  moment  curve  as  the  distance  of  the  adsorbate  from  the 
cluster  is  varied  [17,19].  For  an  ionic  molecule  represented  by 
two  point  charges  +q  and  -q,  with  the  negative  charge  on  the 
positive  z  axis,  the  dipole  moment  is  uz  =  -q  x  r ,  where  r  is  the 
distance  between  the  two  point  charges;  thus  the  slope,  or  first 
derivative,  of  the  curve  is  d:,/dr  =  -q.  For  a  fully  ionic 
molecule  with  q=-l,  dn/dr=-l  and  the  curve  is  a  straigth  line. 
Actually,  in  real  systems,  the  polarization  contributions  can 
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change  significantly  the  value  of  the  slope  because  of  the 

polarisation  of  the  subunits.  This  polarization,  in  particular  of 
the  metal  substrate,  can  lead  to  a  substantial  decrease  of  the 
absolute  magnitude  of  the  dipole  moment.  This  decrease  is  the 
reason  why  changes  in  the  surface  d:pole  or  work  function  are  not 

always  a  good  measure  of  the  ionicity  of  a  chemisorbed  atom  on  a 

surface  [20].  It  is  also  the  reason  why  the  absolute  value  of  t-hi%’.’* 
dipole  moment  is  not  a  measure  of  ionicity  while  the  slope  of  the 
dipole  moment  curve  can  give  a  better  estimate  of  the  bond 

ionicity.  In  general,  the  slope  is  expected  to  be  small  for  a 
covalent  bonding  and  large  for  an  ionic  bonding.  Indeed,  the 
slope  of  the  dipole  moment  curve  as  function  of  the  metal-ligand 
distance  for  CO  chemisorbed  on  Cui  o  and  Cui  4  is  -0.05  and  -0.09 
[15,21],  respectively,  consistently  with  the  covalent  interaction 
indicated  by  the  projection  operator  analysis. 

For  the  systems  here  considered,  Ag4 -?  and  Cuj-O,  the  dipole 
moment  curve  has  been  determined  by  a  Taylor  expansion  about  r=re 
/i(r)  =  Mo  +  Mt  (r-re  )  +  M2  (r-re)^  +  ... 

The  first  term.  Mo  ,  is  )  :  Mi  gives  the  slope  at  r  =  re  ,  while 

small  Mz  values  indicate  linear  behaviour.  For  both  F  on  Ag  and  O 
on  Cu ,  the  analysis  has  been  repeated  by  replacing  the  real 
adsorbate  atoms  with  a  negative  point  charge,  PC=-1  (Table  2). 
The  F  atom  behaves  very  similarly  to  a  -1  point  charge;  both 
dipole  moment  curves  are  linear  and  e.xhibit  similar  slopes  (Table 
2) .  The  slopes  differ  from  being  exactly  -1  because  of  the 
polarization  of  the  subunits,  in  particular  the  Ag4  metal  cluster 
[17] ,  For  O  and  PC  on  Cus  ,  the  dipole  moment  curves  exhibit  large 
negative  slopes  (Table  2),  a  sign  of  dominant  ionic  bonding. 
However,  the  slope  for  the  Cus -O  system  is  larger  than  that 
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obtained  when  we  move  the  PC=-1  along  the  z  direction,  indicating 
that  the  negative  charge  associated  with  the  O  atom  is  larger 
than  -1,  in  agreement  with  the  result  of  the  orbital  projection 
analysis . 

External  electric  field 

A  further  proof  of  ionicity  can  be  obtained  by  applying  an 
uniform  external  electric  field,  F,  normal  to  the  cluster  surface 
[15,21].  In  fact,  an  electric  field  will  produce  large  changes  on 
the  chemisorption  geometry  of  ionic  adsorbates  but  will  have  only 
small  effects  on  covalently  bonded  adsorbates  [15] . 

For  a  cluster  in  a  uniform  electric  field,  the  SCF  first  order 
perturbation  theory  energy,  Ep (F) ,  is  obtained  as  the  difference 
between  the  F=0  SCF  energy  and  _u  x  F,  where  ii.  is  the  field-free 
dipole  moment.  This  is  a  pure  Stark  effect  and  does  not  include 
any  chemical  change  caused  by  the  field.  These  electronic  effects 
are  explicitly  considered  in  the  SCF  variational  energy  in  the 
presence  of  the  field,  EscfCF),  obtained  by  addition  of  the  term 
(I  ri  -  I  Zn  Rn )  X  F  to  the  ' xeld-free  hamiltonian:  i  denotes 
electrons,  N  nuclei,  and  Zn  the  nuclear  charges.  The  fields 
considered  are  +-0.01  a.u.  =  +-5  7  x  10^  V/cm.  The  sign  of  F  is 
such  that  F  <  0  attracts  electrons  from  the  surface  towards  the 
ligand. 

Indeed,  we  found  that  the  field  has  a  large  effect  on  the 
equilibrium  geometry  and  vibrational  frequency  of  F  on  Ag^  .  In 
one  case,  F=-0.01  a.u.,  the  F”  ion  is  pulled  away  from  the 
surface  with  consequent  elongation  of  the  re  value  and  decrease 
of  cue  .  For  F  =  +  0.01  a.u.  the  adsorbate  ligand  is  pushed  against 
the  surface  and  the  metal-ligand  vibrational  frocjuency  increases 
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because  F'  stretches  against  the  ''well"  represented  by  the 
surface  charge  distribution  (Table  3) . 

In  agreement  with  the  picture  of  covalent  bending,  the  field 
has  almost  no  influence  on  the  adsorption  ^eemetry  of  10  adsorbed 
on  a  Cui  4  cluster  [21].  In  fact,  the  equilibrium  bond  distance 

and  the  vibrational  frequency  are  virtually  unaffected  by  the 

* 

external  field  (Table  3). 

For  both  F/Ag  and  CO/Cu,  the  Stark  and  full  S2F  values  of  r* 
and  cue  are  extremely  close,  indicating  that  elctrostatic  erfects 
are  much  more  important  than  chemical  effects,  at  least  for 
fields  of  magnitude  +-  0.01  a.u. 

3.  Consequences  of  ionic  adsorbates  on  m.etal  surfaces 
3.1  Thiocyanate  on  silver  electrodes 

In  this  section  we  summarize  the  results  rf  calculations 
designed  to  provide  a  logical  explanation  of  the  ce.na viour  of  !•'- 
and  S-bound  thiocyanate  species  in  Ag  electrccheni cal  ceils.  To 
this  end  we  consider  a  simple  model.  The  thiocyanate  ton  is 
assumed  to  bind  N-  or  S-down  on  an  on-top  Ag'l^C'  site.  The 
simplest  systems  considered  are  [Ag-NCS]^'"‘  and  [.ig-SCN]  °  ‘  ; 

the  total  charge  of  the  clusters  has  been  changed  to  obtain  a 
very  simple  model  of  the  change  in  electrode  cnarging  as  the 

potential  is  made  positive,  Ag'^ -L- - >  Ag* -b*  .  Ihe  interaction 

of  the  SCI!"  ligand  on  the  top  site  of  an  Ac;  cluster  containing 
one  Ag  atom  in  the  first  and  four  Ag  atoms  in  the  second  layer 
has  also  been  considered  (Fig.  1).  The  m.etal  atc-s  are  fixed  at 
their  bulk  positions  and  the  mietal-ligand  geometry  has  been 
opti.mized  for  II-  and  S-down  coordination  -ede-s.  The  sam.e  1? 
electron  EC?  and  basis  set  [10]  used  for  the  study  of  tne  -E 
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cluster  has  been  eirployed.  In  Ags  ,  the  four  Ag  atoms  in  the 
second  layer  have  been  treated  with  an  ECP  operator  which 
includes  in  the  valence  only  the  5s  electron  [10c] .  C,  N,  and  S 
atoms  have  been  treated  at  all  electron  level;  the  basis  set  for 
C  and  N  is  a  4s3p  contracted  GTO  basis  [18a] ,  the  basis  set  for  S 
is  contracted  to  _^5s4p  [18b]. 

Adsorption  geometry 

For  the  calculation  of  the  molecular  geometry  a  two- 
dimensional  PES  has  been  determined  by  varying  simultaneously  the 
metal-ligand  (M-L)  distance,  r,  and  the  M-N-C  or  M-S-C  internal 
angles,  a.  The  C-N  and  C-S  bond  distanci  have  been  fixed  at  2.20 
and  3.20  a.u.,  respectively.  The  re  and  oe  values  have  been 
determined  from  a  two-dimensional  polynomial  fit.  When  the 
thiocyanate  molecule  is  bound  at  nitrogen  a  perpendicular 
orientation  is  preferred  for  both  Ag-L  and  Ags -L  clusters  (Table 
4).  When  sulphur  forms  the  thiocyanate-surface  bond,  the  molecule 
assumes  a  bent  conformation  forming  a  Ag-S-C  angle  of  about  100® 
(Table  4  and  Fig.  1) .  The  same  bonding  modes  have  been  observed 
in  organometallic  complexes  [22] .  The  nature  of  the  Ag-SCN 
potential  surface  is  very  different  for  the  bending  of  negatively 


charged 

and 

neutral  Ag-L  systems. 

In  particular, 

,  it  is  very 

shallow 

for 

the 

Ag-L-  case  while 

for  neutral 

Ag-L  the  PES 

exhibits 

a 

more 

pronounced  minimum 

(Fig.  2), 

indicating  a 

stronger  directionality  of  the  bond. 

The  N-down  configuration  is  slightly  preferred  in  both  neutral 
and  charged  systems.  However,  the  energy  difference  is  always 
smaller  than  0.2  eV  (Table  4);  this  result  hold  for  both  Ag-L  and 
Ag-l*  clusters.  It  suggests  that  the  two  coordination  modes  can 
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coexist  at  least  in  the  limit  of  low  coverage. 

Bonding  nature 

To  characterize  the  nature  of  the  Ag-thiocyanate  bonding  we 
applied  the  projection  operator  technique  described  in  Sect.  2. 
First  we  projected  the  SCN”  orbitals  on  the  [Ag-NCS]«--i  and 
[Ag-SCN]0'-‘  wavef unctions  for  a  geometry  close  to  equilibrium. 
It  turns  out  that  in  charged  Ag-L  complexes  the  thiodyanate 
ligand  carries  a  full  negative  charge  for  both  bonding  modes  (N 
or  S  down).  For  neutral  Ag-L,  there  is  a  deviation  from  an 
ionicity  -1,  indicating  that  there  is  a  small  covalent 
contribution  (Table  5)  for  Ag-L.  However,  Ag-L"  is  a  better 
representation  for  che.misorbed  SC!I  since  other  Ag  atoms  will 
reduce  the  local  positive  charge  on  the  chemisorption  site  Ag 
atom. 

A  measure  of  the  relative  importance  of  various  charge 
rearrangements  contributions  to  the  Ag-SCn  bonding  can  be 
obtained  by  performing  a  Constrained  Space  Orbital  Variation 
(CSOV)  analysis  [23].  In  particular,  this  makes  it  possible  to 
distinguish  between  intra-unit  pclaritation ,  inter-unit  charge 
transfer  and  electrostatic  effects. 

We  start,  CSOV  step  0,  with  the  superposed  Ag  (Ag* )  and  SCN" 
charge  distributions;  this  gives  the  frozen-orbital  (FO) 
interaction.  Then  we  let  vary  the  Ag  Ag* )  orbitals  in  response 
to  the  presence  of  the  SCri"  anion  (CSOV  step  1)  ;  this  allows  the 
Ag  orbitals  to  polarize  in  the  presence  of  SCN" .  The  next  step, 
(CSOV  step  2)  is  the  Ai  orbitals  variation  in  the  full  space, 
indicating  the  amount  of  charge  transfer  from  Ag  to  SCH"  .  In  step 
3  the  Ag  orbitals  are  fixed  as  they  are  in  step  2,  and  the  SCTJ" 
orbitals  are  varied  in  their  own  space.  This  permits  the  SCM" 

1  a 


unit  to  polarize  in  the  presence  of  Ag .  In  CSOV  step  4  this 
variation  occurs  in  the  full  virtual  space,  so  that  both  SCN"  to 
Ag  donation  or  dative  covalent  bonding  as  well  as  SCN* 
polarization  are  possible.  The  result  of  step  4  is  compared  to 
the  full  unconstrained  SCF  result  and,  if  they  are  nearly  the 
same,  this  means  that  all  important  bonding  contributions  have 
been  considered.  For  each  step,  we  consider  the  interaction 
energy,  Fiat,  with  respect  to  separated  Ag  {Ag*  )  and  SCN~ 
fragments  (Table  6). 

The  nature  of  the  silver-thiocyanate  interaction  in  charged 
and  neutral  Ag-L  complexes  is  basically  different.  In  the  first 
case  the  bonding  originates  mainly  from  the  polarization  of  the 
metal  (CSOV  step  1)  in  response  to  the  presence  of  the  SCN- 
ligand;  of  course,  the  larger  Aga  cluster  is  much  more 
polarizable  than  a  single  Ag  atom  and  the  bonding  in  [Aga-L]", 
1. 2-1.4  eV  depending  on  the  SCN*  orientation,  is  considerably 
stronger  than  in  [Ag-L]"  (0.5-0. 6  eV,  Table  4).  The  Ag 
polarization  dominates  the  bonding  in  the  ionic  complexes ,  but 
small  bonding  contributions  come  also  from  thiocyanate 
polarization  (CSOV  step  3)  and  from  SCN*  to  Ag  charge  transfer 
(CSOV  step  4)  (Table  6). 

For  neutral  Ag-L  complexes  we  start  the  CSOV  analysis,  step  0, 
with  the  superposed  Ag*  and  SCN"  charge  distributions.  This  is 
justified  by  the  previous  projection  operator  analysis  which  has 
shown  the  large  ionicity  of  the  bonding.  Here  the  bonding  is 
dominated  by  the  attractive,  FO ,  electrostatic  interaction,  CSOV 
step  0  (Table  6).  Of  course,  the  polarization  of  the  Ag‘ 
fragment,  step  1,  is  negligible,  but  this  is  an  artifact  of 
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representing  the  r.etal  surface  through  a  single  r.etal  atom. 
Again,  this  contribution  is  expected  to  be  much  larger  for  a  real 
metal  surface.  A  considerable  contribution  to  Ei a t  in  Ag-NCS  and 
Ag-SCN  cames  also  from  the  SCN  polarization  and  from  SCTI  to  metal 
charge  transfer,  CSOV  steps  3  and  4  {Table  5) .  These  effects, 
together  with  the  result  of  the  orbital  projection  iTab_le  5), 
suggest  that  so.me  degree  of  covalency  is  present  -in  this 
essentially  ionic  bend. 

The  results  of  the  CSOV  analysis  help  to  understand  the 
different  bonding  gec.metry  of  the  N-  or  S-bound  surface 
complexes.  We  performed  the  CSOV  analysis  for  various  distances 
and  angles  of  the  [Ag-'ICS]"  and  [.ig-SC!J]'  systems.  keeping  in 
mind  that  the  most  important  contribution  to  Ei  =  are  the  ?0 
Pauli  repulsion,  CSOV  step  0,  and  the  metal  polarication,  CSOV 
step  1,  the  results  can  be  sum.-tarized  as  follows.  In  general,  as 
the  Ag-L  distance  decreases,  both  Pauli  repulsion  and  m.etal 
polarization  increase.  However,  the  Pauli  repulsion  increases 
more  rapidily  than  the  polarization  as  r  becomes  smaller.  Eor  Ac- 
NCS*  a  compromise  between  the  two  tendencies  is  fou.nd  for  r  =  4.5 
a.u.  where  the  potential  energy  curve  e.xhibits  a  mi.nimium.  Per 
a  linear  orientation  of  the  Ag-SCri*  system  we  observe  a  sim.ilar 
competition  between  attractive  and  repulsive  contributions  and 
the  minimum  of  the  potential  energy  curve  is  found  at  r=5.3  a.u. 
because  of  the  larger  spatial  extension  of  the  sulphur  with 
respect  to  nitrogen  charge  distribution.  However,  bending  the 
SCN"  ligand  towards  the  Ag  atom  stabilzes  the  ccmple.x,  while  in 
the  N-down  case  where  this  motion  is  unfavoured.  In  both  cases 
the  bending  moves  the  center  of  charge  cf  the  SC'!'  ligand  closer 
to  the  surface,  thus  increasing  the  metal  polarization ;  however, 
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in  Ag-NCS  this  effect  is  overcompensated  by  the  increasing 
repulsive  interaction  between  the  SCN"  1:t  MO,  mainly  localized  on 
C  and  N  centers,  and  the  Ag  metal  orbitals.  The  overall  result  is 
a  destabilization  of  the  bent  Ag-NCS  structure.  When  SCN  is  bound 
at  sulphur,  the  bending  is  not  accompanied  by  a  similar  increase 
of  the  Pauli  repulsion  because  the  charge  distribution  around  the 
S  atom  is  more  spherical  and  the  repulsive  interaction  between 
the  S  and  Ag  charge  distributions  is  practically  independent  on 
Ag-S-C  angle.  The  result  is  a  minimum  in  the  PES  for  a  =  100°  . 
Field  effects  on  the  thiocyanate  adsorption  geometry 

The  different  geometry  of  N-  and  S-bound  thiocyanate  ligand  is 
the  reason  of  a  different  behaviour  of  the  two  conformations  as 
the  metal  electrode  potential  is  changed.  This  effect  has  been 
simulated  by  applying  an  external  electric  field  parallel  to  the 
Ag-SCN  axis  of  magnitude  +-0.01  a.u. 

The  main  effect  of  the  field  for  N-bound  thiocyanate  is  to 
move  the  ion  up  and  down  depending  on  the  sign  of  F  without 
changing  its  orientation  (Table  7).  When  the  bent  [Ag-SC.N]* 
system  is  considered,  a  field  F=+0.01  a.u.,  such  that  electrons 
are  attracted  toward  the  metal,  further  bends  the  molecule  to 
form  an  Ag-S-C  angle  of  82  degrees.  On  a  real  metal  surface  this 
bending  motion  is  limited  by  the  presence  of  other  Ag  atoms  and 
the  SCN'  ligand  will  lie  flat  on  the  surface.  The  opposite  field, 
a.u.,  orients  the  SCN'  ligand  parallel  to  the  field 
direction  changing  the  shape  of  [Ag-SCN]'  from  bent  to  linear, 
with  consequent  destabilization  of  the  complex,  and  pushes  the 
SCN'  ion  away  from  the  metal  to  such  an  extent  that  no  minimum  is 
found  for  a  field  of  this  .magnitude.  The  same  behaviour  was  found 
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for  the  T.cre  strcr.lgy  hound  [Ag;-SCN]-  cluster;  here  too  the  SCri" 
ion  is  unstable  in  a  field  r=-0.31  a.u. 

Thus.  the  present  results  provide  evidence  that  thiocyanate 
adsorbs  on  Ag  electrodes  at  both  S  and  N  sites  in  the  double 
layer  region;  however,  as  the  potential  is  changed  the  S-down 
configuration  orients  increasingly  towards  the  normal,  and, 
eventually,  dissociates  for  large  negative  fields.  When  positive 
fields  are  applied  the  S-bound  configuration  is  adsorbed  flat  and 
cannot  be  detected  by  IR  spectroscopy. 

3.2  Interaction  of  CO  with  ionic  coadsorbates 

One  of  the  .most  challenging  areas  in  catalysis  is  the 
interaction  occuring  between  an  ad-molecule  and  co-adsorbates 
which  increase  (promoters)  or  reduce  (inhibitors)  the  catalytic 
activity  of  the  metal  substrate. 

Several  theories  [24,25]  have  been  proposed  in  terms  of  simple 
concepts  to  rationalize  the  role  of  promoters,  and  in  particular 
to  e.xplain  the  observed  large  shifts  in  vibrational  frequency  of 
molecules  like  CO  adsorbed  in  the  presence  of  electropositive 
(e.g.  alkali  metals)  cr  electronegative  (S,  P,  Cl,  etc.)  atoms. 
The  proposed  mechanisms  are  essentially  based  on  electronic 
factors  or  on  steric  site  blocking.  Here,  we  consider  an 
additional  mechanism,  purely  electrostatic,  which  is  shown  to 
make  an  important  contribution  to  the  CO  vibrational  shifts. 

The  basic  idea  is  that  an  electropositive  (electronegative) 
adatom  form.s  a  bond  with  the  surface  with  large  ionic  character. 
In  image  charge  theory,  the  electrons  of  a  metal  substrate 
redistribute  in  t.he  prese.nce  of  a  positive  (negative)  charge 
above  the  surface  in  such  a  way  that  the  electric  field  within 
the  metal  is  zero.  This  electron  redistribution  can  be 


represented  by  an  image  charge  below  the  surface  of  opposite  sign 
to  the  "test’*  charge.  The  image  and  test  charges  are  equidistant 
from  the  surface  plane  and  generate  an  electric  field  which  will 
interact  with  the  coadsorbed  molecule.  Thus,  we  model  the  ionic 
adatoms  with  a  pair  of  point  charges.  Specifically,  we  consider 
the  interaction  bf  positive  and  negative  point  charges  (PC)  pairs 
which  are  simple  models  of  ionic  adatoms  and  their  image  charges 
with  the  Cu3 -CO  cluster.  The  geometry  of  the  cluster  with  PCs  is 
shown  in  Fig.  3;  further  details  can  be  found  in  ref.  16. 

Since  our  concern  is  primarily  directed  to  the  C-0  stretch 
vibrational  frequency,  we  consider  geometry  variations  along  an 
internal  coordinate  for  this  stretch.  As  discussed  in  Sect.  2, 
both  Stark  and  full  SCF  shifts  can  be  determined  (Table  3) .  Two 
PC  distributions  have  been  considered;  positive  charges  above  and 
negative  charges  below  the  surface,  hereafter  denoted  as 
(+/-)  PC,  and  the  reverse  case,  (-/+)  PC,  representing, 
respectively,  electropositive  and  electronegative  coadsorbates. 
The  PC  orientation  (+/-)  creates  an  electric  field  which  draws 
electronic  charge  from  the  surface  toward  CO  and  correspond  to  a 
field  F  <0.  This  PC  distribution  induces  large  CO  a>e  shifts, 
about  400  cm"  1  ,  to  lower  frequencies.  The  C-0  distance  is 
substantially  elongated  (Table  8).  These  A<j>  values  represent 
about  70%  of  the  A'ju  e.xperimentally  observed  for  coadsorption  of 
CO  and  K  [26,27].  The  effect  of  the  (-/+)  PC  distribution  is 
opposite  and  raises  the  CO  vibrational  frequency  by  200-300  cm" ‘ , 
thus  reinforcing  the  CO  bonding  and  justifying  a  reduction  of  the 
catalytic  activity.  The  changes  in  re  and  loe  due  to  the  PCs  are 
dominated  by  electrostatic  Stark  effects,  but  they  are 


acccnpanied  by  significant  electronic  rearrangements  mainly  due 
to  a  change  in  Cu  to  2k*  back-donation  [16]. 

The  limitations  of  the  present  model  preclude  our  obtaining 
quantitative  results  for  Are  and  Aioe  ;  moreover,  the  field  created 
by  the  PCs,  about  13  x  IC'  V/cm,  is  probably  too  large  with 
respect  to  a  real  situation.  Nevertheless,  we  have  shown  that 
without  invoking  any  electronic  mechanism,  the  electrostatic 
interaction  alone  causes  remarkable  shifts  in  CO  'i;®  .  This  can  be 
a  first  step  i.n  the  activation  of  the  CO  molecule  accompanied  by 
other,  non  negligible,  electronic  or  steric  mechanisms. 
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Table  1  -  Projection  of  F“  on  Ag^  F ,  -  on  CusO,  CO  on  Cuj-CO, 

for  r=re 


F/Ag^ -F 

O/Cua  - 

0 

CO/Cua 

-CO 

P{2s)  2.00 

P(2s) 

1 . 99 

P(5^) 

1 . 99 

P{2pg)  1.99 

P (2pa) 

1.96 

P  dm) 

4 . 00 

P{2pm)  3.99 

P(2pm) 

3 . 87 

P  (2m*  ) 

0 .32 

AP ( 2pm) 

3 . 46^ 

AP  {2  m*  ) 

0 . 21^ 

Np 

9.98 

9 . 33 

14 . 21 

®  AP((p) 
Cua -L  (L=0 


are  the  values  of 
,  CO)  and  Cua  +  L 


the  difference  of 
subunits  at  r=re  . 


'■P 


projection 


from 


O') 


Table  2  -  Dipole  moment  curves  for  Agi-F,  Ag^-PC,  Cuj -O  and 


Cus -PC  systems  (PC  =  -1  point  charge,  see  text) 


Ag-i  -F 

Ag^  -PC 

Cua  —0 

CU3 -PC 

Mo 

-1.09 

-1.65 

+  0.24 

+  0 . 20 

Ml 

-1.39 

-1 . 32 

-1 . 17 

-0 .72 

M2 

-0.09 

-0.08 

-0.37 

-0.27 
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Table  3  -  Field 

distances  and 

clusters 

induced  shifts  in  metal 

vibrational  frequencies 

-ligand  equi 

for  Ag4  -F 

librium  bond 

and  Cui 4 -CO 

F 

method 

Ag4  ■ 
re  ,  au 

-F 

lUe  ,  cm"  ' 

CUl  4 

re  ,  au 

-coa 

Ij)e  ,  cm"  ‘ 

0 

SCF 

3.31 

277 

3 .95 

184 

-0.01 

Stark 

3.64 

213 

3 . 96 

185 

SCF 

3 .62 

215 

4 . 00 

172 

+  0.01 

Stark 

3.10 

322 

3.95 

186 

SCF 

3.09 

323 

3 . 97 

180 

®  From  Ref.  21. 
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Table  4  -  Optimized  geometries  and  dissociation  energies  of  M- 


and  S-bound 

thiocyanate 

on  Ag  atom  and 

Aga  cluster 

r ( Ag-L) ,  au 

a,  degrees 

De  ,  eV 

[Ag-NCS]~ 

4 . 47 

130 

0 . 64 

[Agj -NCS]- 

4.25 

180 

1.41  ■ 

Ag-NCS 

4.04 

180 

5 . 27 

Ags  -NCS 

4.11 

ISO 

4 . 12 

[Ag-SCN]- 

5 . 39 

107 

0 . 49 

[Aga -SCN]- 

5.09 

105 

1 . 17 

Ag-SCN 

4 . 82 

98 

5.17 

Aga  -  SCN 

4.91 

97 

3.99 

25 


Table  5 

[Ag-SCN] 0 

-  Projection  of  SCN' 

.  -  1  _ 

orbitals 

on  [Ag-NCS]0’- 

■  ‘  and 

[Ag-NCS] - 

Ag-NCS 

[Ag-SCN] - 

Ag-SCN 

Np  (a) 

7.98 

7.93 

3 . 00 

3 . 00 

Np  (:r) 

o 

o 

00 

fl 

7.98 

7 . 97 

7 . 89 

Qp 

-0.98 

-0.91 

-0.97 

-0.89 
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Table  6  -  CSOV  analysis  for  [Ag-NCS]o--‘  and  [  Ag-SC!J]  °  ‘  for  a 

geometry  close  to  equilibrium  (see  Table  4) 


Step 

[Ag-NCS]- 
Ei  n  t  /  iEl  n  t 

Ag-NCS 

El  n  t  /  AEl  n  t 

[Ag-SCN] - 
El  n  t  / AEl  n  t 

Ag-SCN 

Elat/  AEl  a  t 

1 

0  . 

? .  0 . 

-O  .  ?8  .  . 

+3.36/... 

-0.25/... 

+  4.0*2/.  .  .  .  .. 

I 

1 , 

V{Ag;Ag) 

+0.22/+0.90 

+4.09/+0.23 

+0. 24/+0. 49 

+4 . 11/+0. 09 

I 

2. 

V{Ag;all) 

+0.26/+0.04 

+4.14/+0.C5 

+0 . 25/+0 .01 

+4.12/+0.01 

3  . 

V{SCN;SCN) 

+0.42/+0.16 

+4. 82/ +0.68 

+0 . 30/+0 . 05 

+4.41/+0.29 

« 

4  . 

V{SCN;all) 

+0 . 59/+0 . 17 

+5.21/+0.39 

■^0 . 44/  +  0 . 14 

+5 .06/ +0 . 65 

5. 

full  SCF 

+  0.64/ +0 . 0  5 

+5.26/+0.G5 

+0 . 49/+0 . 05 

+5. 16/+0.10 

1 
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Table  7  -  Equilibrium  bond  distances  (a.u.)  and  bond  angles 
(degrees)  for  [Ag-NCS]"  and  [Ag-SCN]"  in  the  presence  of  an 
external  electric  field  F  (in  a.u.).  Only  variational  SC"  values 
are  given 


Field 

[Ag-NCS] 

- 

[Ag- 

■SCN]  - 

re 

a 

re 

V 

0 . 00 

4 . 47 

180 

5.39 

107 

-0 . 01 

4 . 92 

180 

180 

+0.01 

4 . 29 

180 

5.34 

82 

3  This 
F=-0.01 

geometry  is  not  stable  in  a 
a.u.;  the  molecule  becomes 

repulsive 
linear  and 

electric  field 
than  dissociates. 
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Table  8  - 

Equilibrium 

CO  bond  distances, 

re  in  a  .  u  .  , 

and 

vibrational 

frequencies , 

U.-9  ,  in  cm-  ‘  ,  for  CO 

chemisorbed  on 

Cus 

in  the  presence  of  point  charges  (PC) .  Changes  from  the  no  PC 
values,  Ire  and  Au)e  ,  are  also  given. 


re  /  Are 

a)«  /  AuJe 

no  PC 

2.119/.  .  . 

2275/. . . 

(+/-) 

PC 

Stark 

SCF 

2.193/+0.074 
2.206/+0. 087 

1859/-416 

1824/-451 

(-/+) 

PC 

Stark 

SCF 

2.069/-0.050 
2.086/-0 . 033 

2572/+297 

2461/+186 
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FIGURES  CATIONS 


Fig.  1  -  Schematic  representation  of  the  adsorption  geometry  of 
N-  and  S-bound  thiocyanate  ligand  on  the  on  top  site  of  an  Agu 
cluster  model  of  the  Ag  (100)  surface. 


Fig.  2  -  (a)  Potential  energy  surface  for  the  interaction  of  the 
Ag  atom  and  the  SCN"  ligand  as  function  of  Ag-S  distance  and  Ag- 
S-C  bond  angle,  (b)  Potential  energy  surface  for  the  interaction 
of  the  Ag*  ion  and  the  SCN“  ligand  as  function  of  Ag-S  distance 
and  Ag-S-C  bond  angle. 


Fig .  3 

charges 
is  3.70 
Cu  atom 
the  top 
neighbor 


-  Geometry  of  the  Cua -CO  cluster  interacting  with  point 
(PC)  modeling  ionic  coadsorbates.  The  Cu-CO  bond  distance 
a.u..  The  PCs  are  6  a.u.  above  and  below  an  hypotetical 
on  the  Cu(lOO)  surface;  the  distance  of  this  point  from 
Cu  atom  in  Cua -CO  (9.6  a.u.)  is  twice  the  nearest 
distance  in  bull<  Cu. 


